Potential energy threshold for nano-hillock formation by impact 
of slow highly charged ions on a CaF2(lll) surface 

A.S. El-Said, 1 W. Meissl, 1 M.C. Simon, 1 J.R. Crespo Lopez-Urrutia, 2 
C. Lemell, 3 J. Burgdorfer, 3 I.C. Gebeshuber, 1 HP. Winter, 1 J. 
Ullrich, 2 C. Trautmann, 4 M. Toulemonde, 5 and F. Aumayr 6 'Q 

1 Institut fur Allgemeine Physik, Vienna University of Technology, 1040 Vienna, Austria, EU 
2 Max-Planck Institute for Nuclear Physics, 69029 Heidelberg, Germany, EU 
3 Institute for Theoretical Physics, Vienna University 
of Technology, 1040 Vienna, Austria, EU 
4 Gesellschaft fur Schwerionenforschung (GSI), 64291 Darmstadt, Germany, EU 
5 Centre Interdisciplinaire de Recherches Ions Laser (CIRIL), 
Laboratoire commun CEA,CNRS,UCBN,ENSICAEN, 
BP5133, 14 070 Caen Cedex 5, France, EU 
6 Institut fur Allgemeine Physik, Vienna University of Technology, 1040 Vienna, Austria 

(Dated: February 4, 2008) 

Abstract 

We investigate the formation of nano-sized hillocks on the (111) surface of CaF2 single crystals 
by impact of slow highly charged ions. Atomic force microscopy reveals a surprisingly sharp and 
well-defined threshold of potential energy carried into the collision of about 14 keV for hillock 
formation. Estimates of the energy density deposited suggest that the threshold is linked to a 
solid-liquid phase transition ("melting") on the nanoscale. With increasing potential energy, both 
the basal diameter and the height of the hillocks increase. The present results reveal a remarkable 
similarity between the present predominantly potential-energy driven process and track formation 
by the thermal spike of swift (~ GeV) heavy ions. 

PACS numbers: 34.50.Dy, 79.20.Rf, 61.80.Jh 
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The surface topography of many solids experiences drastic modifications when exposed 
to energetic ions. The changes induced depend on the target material as well as on various 
beam parameters such as charge, energy, mass and fluence of the incomingprojectiles and 
can result in well-ordered patterns, such as ripples or self ordered dots . Impact of 

single ions has been demonstrated to induce nano-sized hillocks on metals, semiconductors 
and dielectric targets. Remarkably enough, the hillocks observed have a similar height (a 
few nm) and diameter (20-40 nm) although the beam energies used span several orders of 
magnitude reaching up to GeV. 

Impact of swift heavy ions is known to induce physical, chemical, and structural modi- 
fications not only on the surface but also in the bulk (see e.g. jj, 0, Q and Refs. therein). 
Individual projectiles form cylindrical tracks around their trajectory of a few nanometers in 
diameter. Track formation sets in above a critical value of the energy loss dE/dx of the pro- 
jectiles and occurs particularly in insulators (e.g. polymers, oxides, ionic crystals). Depend- 
ing on the solid, tracks consist of amorphised or defect-rich material. In non-amorphisable 
alkali and alkaline earth halides (e.g. LiF and CaF2) the damage process is governed by 
exciton-mediated defects such as color centers and defect clusters p, Above a critical 
value of dE/dx, damage produced in the core of the track leads to a macroscopic volume 
increase (swelling, 0,[lO]), track etchability Bl, a nd stress jllj]. At the surface of ionic crys- 
tals, swift ions induce nanometric hillocks 12|,[l3j] above a threshold value similar to that for 
swelling 10]. Although numerous experimental data are available for hillock formation due 
to swift heavy ion impact, the principle of the mechanism is still not yet fully understood. 

In this letter we present experiments with slow (v p ~ 0.3 a.u.) highly charged ions (HCI) 
which also induce hillock-like nanostructures on the surface of CaF 2 single crystals. These 
nanostructures closely resemble those created by fast ions. Moreover, we find a strong 
dependence of the formation on the potential energy rather than on the stopping power. 
Most surprisingly, we find a well-defined threshold of potential energy required for the onset 
of nano-hillock formation. Since CaFo is used as an insulator in silicon microelectronic 
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devices 11 1!. 15] epitaxially grown on semiconductor surfaces our findings might be of 
importance for high resolution patterning of thin CaF 2 films on Si and for the creation of 
nanostructured templates for adlayer growth during fabrication of CaF 2 /Si-based epitaxial 
insulator-semiconductor structures. 

Our experiments were performed on air-cleaved CaF2(lll) surfaces. Cleavage is known to 
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result in a fluorine-terminated surface. Contact-mode atomic force microscopy (AFM) in air 
revealed large atomically flat surfaces with occasional cleavage steps separating individual 
terraces. Several freshly cleaved CaF 2 samples were mounted in a vacuum chamber of 
pressure in the 10 -10 mbar range and irradiated normal to the (111) surface with HCI of 
kinetic energy below 5 keV per nucleon. The irradiation was performed at the Heidelberg 
electron beam ion trap Q using «A r - fa = 11, 12, 14, 16, 17, and 18, as wen as »Xe- 
(q = 22, 26, 28, 30, 33, 36, 40, 44, 46, and 48) projectiles during several runs. The extraction 
voltage was 10 kV (for Xe 44+ also 6.4 kV) equivalent to a kinetic energy of 10 kV (6.4 kV) 
times charge q resulting in a projected range between 90 and 140 nm in CaF2, assu ming 
that stopping power and range are unaffected by the high charge state (see below) [18]. 
The beam flux varied between 10 3 and several 10 4 ions/s, and was measured via electron 



emission statistics with close to 100% detection efficiency 19, 20]. After exposure to fluences 
up to (0.5 — 5) x 10 9 ions/cm 2 , the surface of the crystals was inspected in ambient air by 
contact-mode AFM. As reported earlier for CaF 2 single crystals irradiated with swift heavy 
ions, the surface hillocks are stable in atmosphere at room temperature jl^ . 

Fig. H shows examples of AFM topographic images of CaF2(lll) after irradiation with 
Xe 28+ (2.2 keV/amu), Xe 30+ (2.3 keV/amu), Xe 40+ (3.1 keV/amu), and Xe 46+ (3.6 keV/amu) 
ions. Hillock-like nanostructures protruding from the surface are observed for highly charged 
Xe 9+ (q > 30) and fully stripped Ar 18+ ions whereas targets irradiated with X.e q+ (q < 28) 
and Ar q+ (q < 17) projectiles did not exhibit any hillocks. The sharp transition, e.g. between 
q = 17 and 18 of argon cannot be associated with irradiation parameters in an obvious way. 
Moreover, results from measurements with 6.4 • q keV (2.2 keV/amu) Xe 44+ differ by less 
than 5% from the data of 10 • q keV (i.e. 3.4 keV/amu) Xe 44+ ions. It appears that the 
kinetic energy plays no decisive role for the size of the observed nanostructures. 

The AFM images were evaluated with respect to number density, height and width dis- 
tributions of the hillocks. The number of the hillocks per unit area was found to be in 
good agreement with the applied ion fluence, i.e., above the threshold, a large majority of 
projectiles (> 70%) produces an individual hillock. Their height ranges between 0.5 and 1 
nm and their diameter between ~ 20 and 60 nm. Due to the finite curvature radius of the 
AFM tip (nominally 4-5 nm), the diameter (but not the height) of the hillocks is subject to 
a systematic error. The protrusions are rather flat with a diameter to height ratio between 
40 and 60. In contrast to hillocks induced by swift heavy ions ^^|, we observe only a weak 
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FIG. 1: (Color online) Topographic contact-mode AFM images of a CaF2(lll) surface irradiated 
by Wq keV Xe 9+ ions of charge state q = 28, 30, 40, 46. In each frame an area of l/im xi^um 
is displayed. Hillock-like nanostructures protruding from the surface are only observed for Xe 
projectiles with charge state q > 30. Above this threshold, the height and diameter of the hillocks 
increase with ion charge state. 

correlation between the diameter and height value of a given hillock. Furthermore, the size 
data were found to be strongly dependent on the potential energy the projectile carries into 
the HCI-surface collision (Fig. EJ). 

The potential (i.e. internal) energy E p of HCI is equal to the total ionization energy 
required for producing the high charge state from its neutral ground state. E p is known to 
have a strong influence on surface interaction processes such as electron emission, sputtering, 
and secondary ion emission j^j]. For both Xe and Ar ions a remarkably well-defined sharp 
threshold in potential energy (between E p « 12 keV for Xe 28+ and E p « 14.4 keV for Ar 18+ ) 
for hillock formation emerges. Above this threshold, an increase of the potential energy 
leads to an increase of both the basal diameter and the height of the hillocks. Another steep 
increase of the mean hillock diameter potentially indicating a second threshold is found 
between Xe 44+ and Xe 46+ (top of Fig.EJ. 
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FIG. 2: (Color online) Mean diameter (top) and height (bottom) of hillock-like nanostructures as a 
function of the potential energy of Ar q+ (open symbol) and Xe 9 " 1 " (full symbol) projectiles. Hillocks 
are found only above a potential energy threshold of about 14 keV. The error bars correspond to 
the standard deviation of the diameter and height distributions; the solid lines are drawn to guide 
the eye. 



A convenient starting point for an analysis of the observed hillock formation is the inter- 
action of the HCI above the surface involving a series of complex processes on different time 
and energy scales. When the ion approaches the surface, neutralization starts by electron 
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2311 . Deexcitation 



transfer from the target into highly excited states of the projectile 
of the projectile proceeds via Auger-type processes producing primarily low energy electrons. 
Only for very highly charged heavy ions with open K and L shells electron energies up to 
several keV can be expected. For these states, however, radiative decay becomes important 
as a competing deexcitation mechanism with fluorescence yields of typically ~ 12% [24]. An 
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increasing amount of potentia. energy is therefore d 18S1 pated by X-ray emission. Thecritica. 
distance R c from the surface for electron transfer to the HCI can be estimated as 121 



^ y/2qe(8i + e-l) 
Rc (s+l)W ' (1) 

where i is the amount of charge left behind (for the first electron capture % = 1), and W 
and £ are the workfunction and the dielectric constant of the material, respectively. For 
CaF 2 we find R c ~ 0.16^ nm, which sets an upper limit for the time available for the 
above-surface neutralization sequence. As an example, for an ion of q = 40, R c is about 
1 nm and the neutralization time is of the order of 1 fs. As the projectile velocity is also 
proportional to yfq in our experiment, the above-surface interaction time is the same for 
all projectiles with equal acceleration voltage. The transfer of electrons to the projectile 
leaves unbalanced holes in the surface which store part of the potential energy of the HCI. 
It is known from electron-emission yield measurements that ~ 3 q electrons are emitted per 
projectile j^j]. For the impact of a q = 40 ion we therefore estimate a number of about 150 
unbalanced holes (emitted electrons + electrons required for neutralization) created. They 
diffuse only slowly into the material (hole velocity in the valence band derived from tight- 
binding calculations is smaller than 0.33 nm/fs [26]). Furthermore, two holes (F° atoms) 
in adjacent sites recombine to volatile fluorine gas molecules leaving behind a Ca-enriched 
metallic surface. Upon impact of the projectile the target is structurally weakened and 
features fluorine depleted, defect-enriched areas. 

For an analysis and interpretation of our data, we adapt aspects of the inelastic thermal 
spike model developed for swift ions |27[. The underlying assumption is that the initial 
deposition of projectile energy involves the electronic subsystem of the target and proceeds on 
a (sub-) femtosecond scale while the energy transfer to the lattice and the concomitant lattice 
deformation and defect production occurs on a (sub-) picosecond scale. The present case of 
slow HCI differs, however, substantially in two aspects. The primary energy distribution of 
"hot" electrons results from a relaxation process of a hollow atom rather than from a Bethe- 
Born type ionization spectrum. A significant fraction of the potential energy is emitted by 
energetic (~ keV) inner-shell Auger electrons. Moreover, slow HCI deposit their potential 
energy in a shallow surface region, whereas swift ions deposit kinetic energy along the full 
length of their trajectory within a cylindrical volume. 

In the following we estimate the amount of energy and the target volume in which HCI 



deposit their potential energy. From calorimetric measurements it is known |28[ that only 
part of the potential energy is transferred to the target. We suppose that this fraction is 
around 50% with an uncertainty of 20%. The excited target electrons spread their energy 
within ~ 100 fs by diffusion into a hemispherical volume around the impact site with a radius 
large compared to the source volume determined by the hollow-atom relaxation (R c pa 1 nm). 
In turn, the diffusing hot electron distribution transfers energy to the lattice by inelastic 
electron-phonon scattering with a characteristic time constant T e - P h of a few hundred fs. 
Phonon-mediated energy transport leads to further spread and thermalization. 

Assuming, for simpl icity , the same overall thermal diffusion length Ad pa 4 nm as observed 



for swift ions in CaF 2 |27j , the fraction of internal energy Ed is deposited in a hemisphere 
of radius Ad comprising about iV pa 8.5 x 10 3 atoms. If the energy deposition per atom, 
En/N, within this locally heated volume exceeds the melting energy of Em = 0.55 eV/atom 



29] a solid-liquid phase transition is expected. Likewise, for Ep/N > E$ = 1.55 eV/atom 
29j, sublimation should set in. In order to have these energies available at the impact site, 
the HCI needs a potential energy above E\ h j = 14 keV and Eg = 40 keV, respectively. 
Such a crude estimate carries a large error bar of about ±50% due to the uncertainty in the 
effective Ad and the fraction of deposited energy. The estimates are remarkably, maybe even 
fortuitously, close to the observed threshold for hillock formation (Ef}) and for the second 
drastic, almost steplike size increase {E l g). 

It should be noted, that the conceptual difficulty in applying the model of the thermal- 
ization of the internal energy within Ad to slow HCI lies in the fact that the difference in 
internal energy between subthreshold (Ar 17+ ) and above threshold (Ar 18+ ) is emitted in one 
additional K- Auger electron with an energy of Ek ~ 4.5 keV and its large inelastic mean free 
path A^ ^> Ad- The deposition of this energy difference is thus not confined to the critical 
volume of melting or evaporation. Moreover, a large fraction of K-Auger electrons emitted 
near the surface is directly ejected into vacuum and a fraction of K-holes is dissipated by 
X-rays and thus unavailable for thermalization. 

An alternative and additional heating mechanism could be the pre-equilibrium charge- 
state dependent electronic and nuclear stopping in insulators strongly deviating from stan- 
dard values for near- neutral projectiles in equilibrium [30[. Such deviations have been found 
for low but not negligible projectile velocity of v p ~ 0.3 VBohr- In a shallow region at and 
below the surface, a strong enhancement with charge state q of the kinetic energy deposition 
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and, correspondingly, reduction of range is expected for highly charged ions. This could 
increase the energy deposition dE/dx near the surface to above the critical value for phase 
transition observed for swift heavy ions j^. Future experiments at lower v p should shed 
light on the role of this energy deposition process. 

Irrespective of the not yet fully understood details of the heating mechanism, the following 
scenario emerges: hillock formation is the result of local melting and swelling when the energy 
deposition by HCI near the surface exceeds the melting energy Em = 0.55 eV/atom. If the 
energy deposition exceeds the critical value for sublimation E$ = 1.55 eV/atom, evaporation 
should lead to the formation of blisters of enhanced size. Moreover, one should expect the 
transition from blister to crater formation when the evaporation is further enhanced. This 
scenario, however, suggests that crater formation should be more likely for even higher q 
and at near-grazing impact angles when the energy deposition concentrates near the topmost 
atomic layer, and direct evaporation into vacuum becomes possible. 

In conclusion, the bombardment of a CaF2 surface with moderately slow (v p = 0.3 a.u.) 
highly charged Ar and Xe ions produces hillock-like surface nanostructures. The formation 
of these protrusions requires a critical potential energy of 14 keV (Ar 18+ and Xe 30+ ). A 
second threshold characterized by a steep increase of hillock diameter appears at 50 keV 
(Xe 44+ ). In analogy to hillock formation by swift heavy ions, we associate the two thresholds 
with phase transitions of melting and sublimation caused by the deposition of the potential 
energy within the electronic subsystem. The presently discussed scenario suggests future 
investigations of HCI induced nanostructures at smaller v p , larger q, and grazing incidence. 

This work has been supported by Austrian Science Foundation FWF (Projects No. 17449 
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Infrastructure at MPI Heidelberg Germany, supported by Transnational Access granted by 
the European Project HPRI-CT-2005-026015. 
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